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Abstract 

A quantum-mechanical technique is used within the framework of U(2) super- Yang- 
Mills theory to investigate what happens in the process after recombination of two D-p- 
branes at one angle. Two types of initial conditions are considered, one of which with 
p = 4 is a candidate of inflation mechanism. It is observed that the branes' shapes come 
to have three extremes due to localization of tachyon condensation. Furthermore, open 
string pairs connecting the decaying D-branes are shown to be created; most part of the 
released energy is used to create them. It also strongly suggests that creation of closed 
strings happens afterward. Closed strings as gravitational radiation from the D-branes 
are also shown to be created. A few speculations are also given on implications of the 
above phenomena for an inflation model. 
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1 Introduction 



Recombination of D-branes at angles are processes increasing its value recently from 
both phenomenological and theoretical point of view, related to the presence of tachy- 
onic modes which appear as modes of open strings between the two D-branes Ij. In 
the brane inflation scenarios, these are some of the promising candidates to explain 
the mechanism of infiation^j [H] |H] [H] [7j [H] In addition, among various string 

constructions of Standard Model, intersecting brane models are one class of hopeful 
candidates fn\ [ T^ [Tlj [TH] [TB] pTj ^ . in some of which it has been proposed that the recom- 
bination process occurs as Higgs mechanism pTj (see also [HH^n])- On the other 
hand, this system and process can be regarded as a generalized setting of DD system 
and its annihilation process, which has been studied thoroughlv ^^ ^i] P ^j ( ?7] . 

In this way, this is one of the most important phenomena to explore at present in string 
theory. 

The behavior of the system after recombination, however, had been almost unexam- 
ined until recently. Suppose one consider a recombination of two D-branes at one angle, 
wrapping around some cycles in a compact space. One might expect that after the re- 
combination, each brane would begin to take a "shorter cut" and then make a damped 
oscillation around some stable configuration of branes, while radiating RR gauge and 
gravitational waves (and some others), leading to the stable configuration. However, the 
case was that only the shape of the recombined branes at the initial stage was merely 
inferred}^, though the final (infinitely late time) configuration can be determined in 
some cases (e.g. see ref. jHI). 

Recently, in ref. ^Hj, K. Hashimoto and Nagaoka proposed that a T-dual of super 
Yang-Mills theory (SYM) can be an adequate framework to describe the process, and in 
the case of two D-strings at one angle and via classical analysis, they presumed how the 
shapes of the branes develop at the very initial stage* 

Main purpose of this paper is to investigate the process after the recombination via 
SYM "more rigorously" and to reveal what happens in the process. To be concrete, 
we will focus on the two points: time-evolution of the D-branes' shape, and behavior 
or creations of fundamental strings. The former can be of value in that we describes 
time-dependent, i.e. dynamical behavior of curved D-branes though for a rather short 
time-scale (which enables us to discuss e.g. gravitational waves radiated by the branes). 
The former might also be regarded as a first step toward the understanding of branes' 

^For a recent review and other references, see e.g. ref.|18|. 

•I- Just after ref.|2n], the process is also analyzed via tachyon effective field theory in ref.jHIII and very 
recently, the subsequent work was done by K. Hashimoto and W. Taylor in ref.|;-{l|. 
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shapes after the recombination in intersecting brane models as Standard Model, which 
are responsible for physical quantities after the condensation of Higgs fields. The latter 
point is what happens in the process itself. What we mean by "more rigorously" is that 
we set concrete initial conditions and make a quantum-mechanical analysis: These are 
necessary steps in order to understand precise behavior of the system, because the system, 
its tachyon sector in particular, is essentially a collection of inverse harmonic oscillators, 
and its behavior depends crucially on the initial fluctuations of the system, which can 
only be evaluated appropriately via a quantum analysis, but has not been done at least 
in this setting.^ 

The second purpose of this paper is to discuss implications of the above two points 
(phenomena), for the inflation scenario of the setting that two D-4-branes are approaching 
each other at one angle as in ref.|n][Hl, aiming at the understanding of the whole process 
of the scenario. 

For the above two purposes, we consider two types of initial conditions for branes. 
As we will discuss later, in order to make a quantum-mechanical analysis, we need the 
fact that all the modes (except U(l) gauge fields on each of the branes) have positive 
frequency-squareds at t = 0. So, the first one we consider (denoted clS CclSG (I)) is that 
two D-p-branes have been parallel until t = 0, but are put intersected at one angle 9 at 
the instant t = 0. The second one (the case (II)) is that one of parallel two D-p-branes 
were rotated by a small angle 9 and are approaching each other very slowly. The case (I) 
is one of the simplest conditions, rather easier to study the process itself. The case (II) 
is more practical; the case with p = 4 is one of the hopeful setting for the brane inflation 
scenarios jnj[B]- From a practical viewpoint, the case (I) may seem less realistic, but it may 
also be able to be regarded as some local or rough approximation to more complicated 
setting of D-branes. 

The outline of this paper is summarized as follows: 

We consider two D-p-branes for each initial condition and for each p > 1, within the 
framework of a T-dual of SYM, a low energy effective theory of open strings around the 
D-branes when distances between the branes are small, as in ref . [1^1 El] • Detailed set-ups 
are given in section 2. 

In section 2 and 3 we evaluate time-evolution of typical amplitudes of tachyonic fluc- 
tuations in the following way; we consider fluctuations around the D-branes and make 
mode-expansions of them to examine their spectra. Then, defining typical amplitudes 
(absolute values) of the fluctuations in terms of VEV's of their absolute square, we obtain 

§In the case of inflation and other systems, see ref . [32] [23 j a-nd for a non-perturbative analysis, see 
ref.IHU. 
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them by using a time-dependent wave function of each mode and focusing on contribution 
of tachyonic modes. We note that though we use WKB approximation (second order one 
in fluctuations), it is enough good in the cases until the tachyons blow up, because higher 
order terms are suppressed by a factor of 0{gs), as we will discuss later. 

In section 4 we discuss time-evolution of D-brane's shape after recombination; we show 
that a seemingly queer behavior of the shape appears and explain why it happens. The 
information of the shape is obtained from a transverse U(2) scalar field by diagonalizing 
its VEV's (i.e. choosing a certain gauge) and looking at its diagonal elements jSH]-^ We 
denote this as "physical gauge" because in the gauge direct correspondence holds between 
d.o.f. of each open string starting from and ending on certain branes and d.o.f. of each 
U(2) matrix element. 

Then, in section 5 we argue that open string pairs (composed of strings with opposite 
orientations) connecting the two decaying branes are created. Our argument is based on 
the behavior of the U(2) gauge and scalar fields in the physical gauge; investigation of 
its non-vanishing energy density, the correspondence of d.o.f. stated above, and charge 
densities of two U(l) symmetries on the D-branes lead to such conclusion. We also evaluate 
number density of the open strings and distribution of energy to the open strings, and 
show that most part of the released energy are used to create the pairs. 

In section 6, we discuss creation of closed strings, going out of the framework of 
SYM: We show that the decaying branes radiate closed strings as gravitational wave, 
regarding the branes as a collection of massive small objects; its energy flux is evaluated 
quantitatively. We also discuss the system after some relatively long time, i.e. beyond 
the reach of SYM; creation of closed strings are strongly suggested after the two branes 
are separated further. 

In section 7, we discuss implications of the above results for a brane inflation scenario, 
on the efficiency of the reheating and mechanism to create "reheated" gauge fields and 
fermions. 

Finally in section 8 we present summary and discussion; we discuss some points on 
non-commutative scalar fields appearing in the case (II). In the appendix we give a short 
review of giving a width of a Gaussian wave function for each mode at arbitrary time. 

case (II), there appears a non-commutativity between VEV's of scalar fields, which implies their 
uncertainty relation. We will discuss this point in section 8. 
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2 Set-ups and preliminaries 



In this section we present set-ups and preliminaries for a quantum-mechanical analysis 
for the case (I) and (II) . 

First we present our notation and set-ups; in this paper we set Is — 1 for convenience 
(and revive it when needed), and set the string coupling gs very small and consider D- 
p-branes in a 10-dimensional fiat target space with a metric g^^i^ = diag(— 1, 1, .., 1) and 
coordinates xm {M — 0, ..,9). The world-volumes are parametrized by xo, xi, ■ • • Xp. 
We compactify space dimensions Xp, ■ ■ ■ ,Xg on a (lO-p)-torus with periods for a — 
p, • • • , 9. We denote Xo as t and Xp as x below. The initial condition for the case (I) is the 
configuration that two D-p-branes have been parallel to each other with a displacement 
yo until t — but are put intersected at one angle 9 at the instant t — 0. The setting is 
depicted in Fig.l. 



(i) t < (ii) t = 

Xp+i Dp Xp+i Dp 




Lp Dp Lp Dp 



Figure 1: The initial configuration for the case (I) 

The initial condition for the case (II) is the one that one of parallel two D-p-branes were 
rotated by a small angle 9 and are approaching each other with a small relative velocity 
V. The simplest setting for each case is represented in Fig. 2. We note that x is parallel 
to neither of the branes, but can parametrize a world-space of either branes, though the 
energy density is + (3'^ times larger then the parallel embedding (where (3 = tan{9/2)). 
Although there are multiple intersecting points in cases of compact dimensions, in this 
paper we assume that the recombination happens only at the point x = for simplicity. 

Let us give a brief review of the process. The configuration represented as (ii) in Fig.l 
cannot be distinguished from that of the two branes in which both are bent and touched 
at Xp = Xp+i — since they have the same energy density and fiux. The latter is unstable 
and each brane begins to take a shorter cut to save the energy. This is the recombination 
process. Their shapes just after that are expected to be like a hyperbola. This instabihty 
leads to the existence of a complex tachyon field with a mass-squared 

= -9, (2.1) 
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Figure 2: The initial configuration for case (II) (t=0) 



appearing as modes of a string connecting the two D-branes, as shown in ref.J^. We 
will investigate the behavior of the system using T-dualized versions of U(2) SYM as in 
ref.jSE][2n!. 

For each p, the SYM action with Xp+i, ■ ■ ■ ,xg T-dualized takes the form 

S2Dp = -Tp J (F+^x Tr{ (1+) \f^uF^' + - [X„ X/} (2.2) 

where Tp = I/qs is the tension of a D-p-brane. F^y is the field strength of the world- 
volume U(2) gauge field for yU = 0, 1, Xi ioi i = p + 1, .., 9 are U(2) adjoint scalar 
fields corresponding to coordinates transverse to the branes. We note that to apply 
SYM, the relation Lp+i ^ LpQjl < ls{= 1) should be satisfied, leading to the constraint 

Lp < ^. (2.3) 

We discuss dynamical behavior of the system for t > 0. The "background" D-p-branes 
for each case for t > are represented by the configuration (VEV's): 

Xp+i=(^^Q Xp^2 = --- = Xs = A^ = 0. (2.4) 

where P = tan(^/2). For the case (I) {z{t) = 0), this is T-dual to the configuration of 
two D-(p+l)-brane with a constant field strength -f'p.p+i = P, as discussed in ref . PU] . 
Actually, in ref.|35j, tachyonic modes around the T-dualized (-Fp,p+i = P) background 
were shown to appear in off-diagonal ones of Ap and Ap+i in a certain gauge, which is 
T-dual to Xp+i in the present case. However, the connection via T-duality does not of 
course means that the two systems are physically equivalent. (Actually, as we will show 
later, T-duality has a non-trivial effect for the system to change of the width of a wave 
function by a V2.) Thus, we examine the spectrum from the beginning. For case (II) 
{z{t) 7^ 0), z/2 is a VEV of an adjoint Higgs field which decreases slowly. U(2) gauge 
symmetry is spontaneously broken to U(l) x U(l) by the non-trivial VEV's. Since we 
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set Qs and 6 very small, the force between the two branes is so weak[^ j2Hl(and for short 
distance, see 0) that we can regard the velocity as a constant, leading to 

z = zq — vt. (2.5) 

In this case potentially tachyonic modes should appear in the off-diagonal elements of Ap 
and Xp^i again. In addition, there might be a possibility that tachyonic modes appear in 
off-diagonal elements of Aq and Xg because the two D-p-branes approaching with a small 
velocity is T-dual to a two D-(p-|-l)-brane system with a constant field strength (though 
it is electric), and is also considered to be a Wick-rotated version of branes at one angle 
7;|28]. So, we denote the fluctuations, including those of Aq and Xg, as 

^^^C 11 ^--(°o o') (-) 

(and denote in particular, Cp as c), and discuss their spectra below. 
The part of the action second order in the fluctuations is 

S2Dp\2nd = Tpj (F^\[-d^cid^d' + \d^c^\^-d^(rd^d-d^w*d^w 

+2i(5{c*d - cd*) - 2iPx{c*^d>'d - c^d^d*) + iv{clw - cqw*) 
-iz{cld^w - c^d^'w*) - Wxf + 2')|c^|' - \2l3xw - zd\% (2.7) 

We consider such a gauge fixing condition that it is written as Dj^^A^ = df^A^—i[A^^\ A^] = 
in a T-dualized pure Yang-Mills system, which gives 

d^,c^ = i{2(3xd + zw). (2.8) 



This is a generalized version of the condition chosen in ref . joS] jdo] . Integrating by part 
some terms and using ()2.8|) . we have 

S2Dp\2nA = Tpj d^^+^x[~cl{-{d,f + A{Pxf + Z^}C, - d*{-{d^f + A{Pxf + Z^}d 

-w*{-{d^f + A{[3xf + z^}w - Ai(3{c*d - cd*) - 2w{clw - CqW*)]. (2.9) 

With the redefinition of the fields as 

^_c + id ~ d + ic 

the term —4:if3{c*d — cd*) in ()2.9j) give rise to their mass terms — (— 4/3)|cp — {Af3)\d\'^. 

Here we set z{t) = v = and discuss the case (I) for a while. Let us consider 
mode expansions of the fields c, d and Cf^ (/t = 0, 1, ■ ■ ■ ,p — 1) with respect to x and Xp, 
[fl = 1, ■ ■ ■ ,p — 1) (the latter for D-p-branes with p > 2) to examine their spectra. As 
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for x^, we use Uk{xp) = e^''^-^'' , free wave functions with a wave vector kp,. As for x, the 
operator for each field is written as 

[-{d,f + A{(3xf-Aa(3] (2.11) 

where a = 1, —1, for c, d and cp, respectively. If it were not for the periodicity of x, the 
eigenvalues and eigenfunctions could be given by and /„ as 

ml = 2p{2n-a) (2.12) 
fn = -^(^f"e-^-'H^{^Px) (2.13) 

for each non-negative integer n where are Hermite polynomials; since x direction is 
compactified, exact eigenfunctions are expected to be mathematically complicated, e.g. 
related to theta functions as in ref . [221 ISni • However, we can see that /„ in (j2.13|) is 
localized around x = with a width 5 = 1/ \/2p ~ 1/ \/9. So, if one consider Lp such that 

^ « (2.14) 

one can use (j2.12p and (|2.13|) as an approximate set of eigenvalues and eigenfunctions. 
Then, we can see that the lowest mode of c (we denote as Cq) is a complex tachyon field 
with a mass-squared 

(mo)' = -2/3, (2.15) 

which agrees with ()2.H) obtained in ref. for a small 9, while the lowest modes of Cp are 
massless and the others are massive. We note that the condition ()2.14|) is compatible with 
the one ()2.3p if ^ <^ 1. So, we consider such a range of Lp throughout this paper. (We 
set 6' ~ 2/5 = 10~^ and Lp = 100 later.) We also note that the approximate width of /„ 
is larger than the eigenfunction in the T-dualized theory in ref.jnE]|Sn] by a factor 
which may be one of non-trivial effects of the T-duality. 
Expanding c and the other fields such as 

c{xt,) =J2j ^^^^y^c„,fc-(t)Mfc(x^)/„(x), (2.16) 

we obtain the action 

where lj^ ^ for t > is given as 

uj^^} = k^ ^2i5{2n-a) (2.18) 
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where k = y{kfiy. (Note that a = 1,-1,0 for Cn,k,dn,k,Cu,n,k, respectively.) This looks 
much like the action of a collection of harmonic oscillators with a mass Tp, although it 
includes modes with negative frequency- squareds; we can easily see that the modes co,a; 
with < k'^ < 2(3 have negative tu^'s (while all the other modes have non-negative c<j^'s). 
Thus, they are the tachyonic modes in the case (I) with p >2. For the case of D-strings, 
there is no d.o.f. of fc^ and cq is the only (complex) tachyonic mode with (ujq)^ = —2/3 for 
t > 0. 

For t < 0, one should set 2(3 = 6 = and instead turn on a displacement of the branes 
Ho in Xp+i direction, i/q contribute to the second order action as a mass-squared for each 
of the off-diagonal fluctuations. Since the direction of Xp is compactified, the momentum 
component kp is discretized as kp = 2'Knp/Lp for an integer rip, and the mode functions is 

gnM=^'''V^v (2.19) 
The modes c-^^j^ = {c^^^^^ + id'^^lf}) / ^2, corresponding to ^„p(x), have such {io^^lkf as 

^"l ' = {yof + {kf.f + {2'nnp/Lpf (2.20) 

(where we attach the index (0) to the qualities defined for t < 0). 

Let us move on to the case (II) where 2 7^ 0. As seen in ()2.9|1 . merely plays a role 
of a mass-squared for each field. Making mode-expansions of c, ci, and c/i in terms of Uk 
and /„, we have the action ()2.17|1 with frequencies 

^'n,k ^ = k^ + 2/3(2n - a) + (^o - vt^ (2.21) 

where a = 1, —1, for Cn,kdn,k, Cu,n,k, respectively again. As we will discuss later, we need 
the fact that all the modes have positive frequency-squareds at t = 0, so we set zq = 0.5/^, 
relatively a large value but not to exceed the applicable range of SYM. Since z{t) decreases 
gradually to zero, again (potentially) tachyonic modes are the modes Co,fc with k"^ < 2(3 for 
the case with p >2, and Cq is the only (complex) tachyonic mode for the case of D-strings. 

Qualitative behavior of the system in the case (II) is as follows: at t = 0, all the modes 
have positive w^'s, hence have some small fluctuations. As z{t) decreases to z = y/2(3 
(which gives lu'q q^ = 0), Co,o starts to blow up, and so do the modes Cq^^s with a small \k\ 
one after another when z becomes a critical displacement zl which gives lu'q ^ = 0. The 
blow-ups together give rise to the tachyon condensation (to cause recombination). 

In the next section we will focus on the potentially tachyonic modes Co,fc with 2/3 S> v, 
and make a quantum-mechanical analysis using the action (j2.17p with the obtained tu^'s. 

We note that one can regard Cq and w as non-tachyonic fields for the parameter 
2(3 ^ V. As for the fields Cq and w, the last term in (j2.9|) might also seem to produce 
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mass terms of linear combination of Cq and w, in the same way as the case of c and 
d, with 2(3, c,d replaced by v,Co, w, respectively, resulting in a complex tachyonic field. 
However, the sign of the kinetic term of Cq is different from that of w, and, if one carry out 
the diagonalization, it results in a non-local action including square roots of differential 
operators; it is unclear whether cq and w include tachyonic modes. Let us consider here 
the case with 2(3 ^ v. Then, even if there would be tachyon modes composed of cq and w, 
the critical distances z'Js of them are much shorter than those of c, and it is much before 
the postulated tachyons of cq and w get "awake" that the tachyon modes of c condensate 
and the recombination proceeds. Thus, we consider the parameter region 2(3 ^ v and 
deal with cq and w as non-tachyonic fields. 

3 A quantum-mechanical analysis on tachyon con- 
densation 

In this section we analyze in detail how the tachyon fields condense, applying quantum 
mechanics to the action ()2.17|1 . We assume that the system is initially at zero temperature, 
i.e. at the ground state, for simplicity, though the essential point of our argument will 
not be changed as long as the temperature is low. 

One of the two basic ingredients (we will use here) is that in a harmonic oscillator with 
a positive constant uj"^, the ground state wave function is Gaussian, and its variable q has a 
typical amplitude (absolute value) of fluctuation gtypicai = V< > = i/V2muj = A/^/2 
where < A > denotes a VEV of an operator A, m is a mass of q, and A is a width of 
g's wave function. Another basic one is that even if u"^ is time-dependent and becomes 
negative afterward, time-evolution of its wave function and hence that of Atypical can be 
deduced within WKB (second order) approximation, by using its propagator or kernel, 
as in e.g. ref.jHZI- Taking these two into account, time-evolution behavior of the typical 
amplitude of each mode can be evaluated using quantum mechanics, if an initial condition 
is such that u'^ of the mode is positive at t = and change into a negative one later: the 
initial typical value is \J gs/2uj{t = 0) since in this case the tension Tp = 1/gs corresponds 
to the mass. Since the wave function at t = is Gaussian, the absolute value of the wave 
function at any t remains to be Gaussian within WKB (second order) approximation, and 
the system is characterized only by its width A(t), which is equal to a typical value of 
the variable g(t)typicai times ^/2. (These are reviewed in the appendix A. See ref. [SZI for 
more detail.) We note that typical amplitudes of the modes with non-negative a;^'s are 
suppressed by ^/g^ in their numerators as |cn,fc(t)|typicai = \l Qs/'^^if) {n > 1). We also 
note that the WKB approximation is quite good until tachyonic modes blow up because 
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higher order terms are also suppressed by an extra factor of 0{gs). 



If we define a typical amplitude of fields as |c(x^)| typical = V< 1^(3^/^)1^ >, "we have to 
sum \cn,k\ty^icax{.fnY\un,kf' with respcct to n and integrate it with respect to kp_ except for 
the case of D-strings, but the integral is divergent. However, this is a kind of divergence 
that appears in a usual quantum field theory and has nothing to do with the blow-up of 
tachyons. Thus, we ignore the contribution of all the non-tachyonic modes and define the 
typical amplitude of c(x^) in the case of D-p-branes with p > 2 as 

|5(x,)kypical = { l^^^^ J-^, < |5o,.(t)r > {hfY^' = (3.1) 

where A{t) is a time-dependent part of \c{x)\ and /o(a;) = (2/5/7r)^/^e~^'^^. In the case of 
D-strings we define 



|£(a;^)|typicai = V< |5oWP > = A{t)U (3.2) 

In this paper we regard the quantity p.ip or ()3.2p as a VEV of |c(x)| due to tachyon 
condensation. As in ref.[37J and the appendix A, the time-dependent width Afc(t) of the 
wave function of Co,fc(t), which is equal to -\/2|co,fc(^)|typicah can be written in terms of the 



initial width Afc(O) = y gs/u!k{t = 0) and two independent solutions to the equations of 
motion 

= -^o,kCo,k- (3.3) 

We denote the solutions as "basis functions" . Below we solve the equation of motion for 
each case and use the basis functions and A^i^O) to evaluate |c(x^) | typical- 

For the case (I) with p > 2, we can easily see that the basis functions for t > are 
cosh(i^t) and sinh(i^t) where K = y^— cJq ^ = \/2f3 — k'^ (> for tachyonic modes), and 

the initial width for {t < 0) is Ai°),(0) = v^l^i^^t = 0)|typicai = where 
ool^^^k is given in (j2.2U|) . However, the mode (eigen-) functions (of the operator with respect 
to X arising in the action) change at t = from gnp{x) to fn{x), so we have to "translate" 
l^o°fc(^ ~ 0)|typicai into \co^k{t = 0)|typicai- (We note that Co,a: is related only to c^^J^^ with the 
same since the system keeps translation symmetry in directions even at t = 0.) If 
we define N^^np = I dxf*gnp, it holds g^p = En N^^npfn and we have Cn,k = Hup Nn,npCnl,ki 
and hence 

Afc(t = 0)2 = 2\coAt = 0)|t'ypical = E \No,np\'^. (3.4) 

rir, UJ„ u 



With the Afc(t = 0)^, the typical value is given through ()A.6|1 as 

2A,(0)2ir2 



co.fcltypicai = = — 7. — cosh Kt + smh Kt, (3.5) 
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and we can obtain the typical amplitude |c(x^)| as 

i£(a^.)i?yp.ai = {A{t)for = jf' k^-'dk^{Mx)r (3.6) 

where flp-2 is the volume of (p-2)-sphere. This is the formula for time-evolution of the 
typical amplitude of the tachyonic fluctuation for the case (I) with p > 2. For p = 1 
(D-strings) case, we have 

|5o(x,)LV,i = {A{t)for = ^^^(/o)^ (3.7) 

where Ak{t) is given in ()3.5|) . 

For all value of p, order of \c{x^) |typicai in terms of gs is 0{gs). In addition, \c{x^) |typicai 
has x-dependence through the Gaussian function /o(a;). This means that the excitation of 
the tachyon is localized around a; = as in ref . jSSl EH] |2H] though the width is relatively 
broad such as 5 = 1/ ^/J3. 

A rough behavior of A{t) can be read from the factor eV^^^^~^* ~ gv^* in ()H.fi|l or 
flH.5|l . This can be inferred from the classical solution to the equation of motion, and a 
typical time scale is defined as T^^^ = l/y/2j3. The formula for A{t) is a bit complicated, 
but if we set Lp = lOOZ^, we have J2np=-2 l^n.npP — 0.9987, so it suffices to have the 
terms with Hp = 0, ±1, ±2 in ()3.4j) . If we set yo = 0.5/^, relatively a large value, we can 
approximate uJnp,k — Vo for the n^'s, and hence Afc(t = 0)^ ~ Qs/Uo- So, we have 

|5o..L'ypicai ^ ^i^pZT^ si^h^ - kH + ^ cosh^ pp-kH}. (3.8) 

Since uq 3> 2(3 now, the first term in ()H.8|1 is dominant. The time-dependence of y4(t)'s 
for p = 1,2,4 is plotted in Fig.3, for a specific value of the angle 6' ~ 2/3 = 0.01 and 
Qs = 0.01, which we set on the basis of the brane phenomenology (e.g. refs.j2]|H])- We 
can read a tendency from ()3.6p or Fig.3 that the blow-up is more delayed for the case 
with a larger value of p. This fact can be explained as follows: the tachyonic modes with 
less momentum tend to blow up more rapidly, but there are larger number of tachyonic 
modes for larger momentum. That is, except for the case of D-strings, for larger value of 
p, the most "active" tachyonic mode which has = is more smeared in some region of 
k and the start of the blow up becomes slower. 

Next we derive A{t) for the case (II). The frequency-squared for co,fc is 

tOo,kity = izo-vtf + kl-2p. (3.9) 

We note that each cUqI is positive at t = for zq = 0.5/s. The behavior of the system in 
terms of its wave function is as follows: As z{t) decreases gradually, so does each uJq 1 and 
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Figure 3: logi^ A(t) vs. t/T*^^).The topmost graph is for D-strings, the middle is for D-2- 
branes and the lowest is for D-4-branes, with ^ ~ 2/3 = 0.01, gs = 0.01 and T^^^ = lOlg. 

the wave function co,a,.) of co,a; spreads little by little. At t = to = {zq — \^2P)/v, cUqq 
(with the lowest momentum k = 0) becomes zero and then negative, i.e. tachyonic, and 
the corresponding wave function begins to spread radically, resulting in the blow-up of 
the mode. After that, the modes with lower k are also getting tachyonic and their wave 
functions also start to spread one after another, leading to the tachyon condensation. 

To obtain A{t) in this case, we have only to find two solutions to ()3.3|1 with ()3.9|1 and 
substitute it for ()A.6|) . However, its exact solutions is a bit difficult to deal with, so we 
find some approximate solutions: Let us define t^{k) = {zq — y/2j3 — k'^)/v which gives 
^^o,k{t^) = 0, and rewrite the cJq^ using t^: as 



where we denote the coefficient = \/8v'^(3' and [3' = (3 — k"^ /2. If t approaches t^, the 
second term quadratic int — t^ of ()3.10p becomes so small as to be neglected compared to 
the first term. (We denote the time as t = ti.) Thus, for t > ti, two Airy functions can be 
used as basis functions. On the other hand, taking a closer look will let one notice that 
(13. 3p with (j3.10|) has the same form as the Schrodinger equation in one dimension with 
relatively slowly- changing potential (if we replace thjx). So, when ujQ^k{t) ^ is relatively 
large and changing slowly, we can apply "WKB approximation" to get two approximate 
basis functions. Thus, if this approximation can be applied until the time t = ti, we 
can follow the time-evolution of \I^(t, Co,fc) and Afc(t) by using the "WKB-approximated" 
solutions until t = ti, and then letting the Airy functions take over the role of basis 
functions. We take this prescription to evaluate their time-evolution. 
The "WKB-approximated" basis functions are 





(3.10) 




tl 
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02(t) = ^o~fc^' sm( / dt'u;o,k{t')). (3.11) 
This approximation is valid as far as 



^ < ko,fcP and l^^^l -C |u;o,feP- (3.12) 
On the other hand, the condition for the Airy functions to be used is 

. (3.13) 

V 

So, if we choose ti = t^: — y/SP'/lOv to satisfy ()3.13p (please note the sign), the condition 
to satisfy ()3.12|) is written as 

V < 2/3'. (3.14) 

which is compatible with the condition to regard cq and w as only massive fields. So, in 
this paper we consider such a system. One may worry that tachyonic modes with large k 
do not satisfy ()3.14|1 for any set of /5 and v. In fact, in a set-up with v <^ 2/3, v is so slow 
that the tachyonic fluctuation c blows up enough much before the tachyon modes with 
large k "wake up" . So, the above point is not practically a problem. 

Using the Airy functions Ai{T) and Bi{T), we obtain the width Afc(t)^ of |\['(co,fc, t)P 

as 

= gA-{B^'{n)A^{T)-A^'{n)B^{r)y 



UJi 



2 



+ ^^^^^{5^(ri)A.(r) - M{n)Bz{r)r], (3.15) 

where Ai'(r) = r = al^^{t - t*), Ui = uJo^k{t = ti){= and ti = al^^{ti - 

t^,). And we can obtain |c(x'') |typicai for p > 2 for case (11) by substituting ()3.15|) for ()3.6|) . 
For the case of D-strings, |c(a;'')| = (Afc=o(^)/o)^/2 where A^^oit) is given in ()3.15p . 

As in the case (1), the blow-up of the tachyonic fluctuation is localized around a; = 
with the width S = 1 14^. A rough behavior of A{t) is ~ exp[|{(8v2/?)i/6(^^ _ 1^)^% so 
we can see that a typical time scale of the decay is T^^^^ = l/{8v^py/^ = l/{4v^9). That 
is, the typical time scale depends more strongly on the relative velocity of the branes 
than on the angle 6. The reason for this is that after reaching the critical distance, the 
larger the velocity is, the earlier potentially tachyonic modes "wake up" to become truly 
tachyonic and blow up. Time-dependent behavior of A{t) in the speciflc case of D-4-branes 
with 2(3 = IQ-^ and v = 10"^ is plotted in Fig.4 where T^^^^ ~ 36.84/,. 

In this way, we have obtained, for both cases, explicit functions of typical amplitudes 
of tachyonic fluctuations |c(x^)|typicai = ^(^)/o- 
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Figure 4: logio ^(t) vs. {t-to)/T^^^^ for D-4-branes; 2/3 = 10"^ v = 10"^ and gs = 0.01. 

For convenience, let us define the following quantity for a function of k, g{k): 

« 9{k) y J^,^9ik) (3.16) 

where it is normalized as << 1 >>= For a later use, we also define a "typical 

amplitude of tachyon's wave vector" of tachyonic fluctuations as = (<< k'^ » 
I << 1 >>)^/^. By v^P we can estimate which momentum-possessing tachyon modes 
contribute the most to the blow-up, since the weight for a tachyon mode with a momentum 
k in is (Afe)V2, as in dSU with |co,fc|typicai = Afc(t)V2. 

4 Time-evolution of the D-p-branes' shapes 

In this section we discuss time-evolution of the D-p-branes' shapes and show seemingly 
queer behavior of the brane's shape occurs. Then, we give a physical interpretation of the 
behavior to conclude that it is not actually queer. 

With ()2.10|) and ()2.16|1 . the typical value of d{Xfj^) is 

rf(a;^)typical = -z^fel|P^ = -1 ^^^^^ A_^Co,kit)typic^lUkiXp)foix) (4.1) 

since ^typical is negligible. We note that each Co,fc(t)typicai has its phase as Co,fe(t)typicai = 
|co,fc(^)|typicaie*''''', and so does (i(x^)typicai- We regard this quality as the VEV of the field 
d{x^). Then, the transverse scalar U(2) field Xp+i has the VEV of the form: 

- ( C '^r ) ■ '^•^> 

Let us get geometrical information on the branes from ()4.2|1 . first for the case (I). In 
this case there is only one transverse scaler field which has non-trivial VEV elements, and. 



15 



all the scalar fields are commutative with each other. Thus, as done in such conventional 
cases and in ref. it is a logical step to choose the gauge to diagonalize the VEV's of 
the scalar field and interpret the diagonal parts as the positions of the two branes. Xp+i 
is diagonalized as 



Xp^i — UoXp+iUq ^ — ^ 



y{x,t) 
-y{x,t) 



(4.3) 



where y = ^ {(3xy + |(itypicaiP- Uq is the transformation matrix to diagonalize Xp+i, whose 
explicit form we present here for a later use, as 



TT ^ 1 / + /3x|d |typical ^typical ) (A A\ 

^0 V2^\d\„ y _y^Zr^4yp.^^i y^T^MItypical )' ^ ' ' 

Since it holds |(itypicaiP = |c(a;^)typicaiP/2 = |c(a;^)|typicai/2* which we have obtained in the 
previous section, the formula for the shape of one of the recombined branes is 



y{x,t) 



\ 



(/3x)2 + ^^W^e-2/3-^ (4.5) 




This is the same form as obtained for D-strings in ref. |29]. However, we can make a 
precise analysis on the shape of the branes at arbitrary t based on ()4.5|) . since we have the 
explicit function of A{t) for each case including the overall factor and fine coefficients; We 
can know explicitly when and which applicable condition of the approximations (SYM, 
WKB and some others) breaks. We will investigate the behavior of the system in detail 
using ()4.5|) with A{t)^s obtained in the previous section. 

For the case (II), its essential difference from the case (I) is that there are two scalar 
fields which have non-trivial VEV's of Xp+i and Xg. Since the two do not commute each 
other after tachyon condensation, it is probable that one cannot determine definite posi- 
tions of the branes in both of the dimensions Xp+i and xg simultaneously (i.e. uncertainty 
relation). Here we assume that we can interpret each of Xp+i's diagonal elements as the 
position of the two branes in the Xp+i direction if we focus only on the positions in Xp+i 
direction regardless of x'^^ direction. We will return to this problem in the final section. 

At the initial stage, the shape of each brane is approximately one part of hyperbola 

as in Fig. 5, which is consistent with our intuition; though typical is dependent on x 

through the Gaussian fo{x), its width l/y/P is broad and its differential coefficient is so 

small. The larger |(i|typicai becomes, the further the distance of the two D-brane becomes. 

*In the last equality we use the approximation |cQj,(t)typicaiCofe' (i)typicaiw^Ufc' | = {2t:)p^^S^'p^^'> {k — 
k'){Ak)^ /2 where Uk — e^'^'^^'^ and A^: is given in (|3.5|l . This holds good because only the contributions 
of Cfc and Ck' with k = k' has a coherent phase factor while the others do not. 
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Figure 5: The shape of D-4-branes; 61 ~ 2/3 = 0.01 and A{t) = 0.01 



That is, I (i I typical corresponds directly to such a relative motion between the recombined 
branes as discussed in ref . [21] . 

So far, the shape of the p-branes is within the reach of our intuition. After more time 
passes, however, a seemingly queer behavior occurs; the shape of each brane deviates 
from the approximate hyperbola and comes to have three extremes due to localization 
of tachyon condensation. This happens for the cases of D-p-brane for p > 2, and is also 
expected to happen for the case of D-strings, as we will show below. ^ 

The values of x giving extremes of y{x,t) are formally given by 



which means that each brane's shape has three extremes if it holds 

A > (7r/5/2)i/^ = A,,iticai. (4.7) 

Thus, all we have to do is to show that the condition (j4.7|) is compatible with applicable 
conditions of the approximations we have used. One of the important conditions is that 
of the WKB approximation that higher order terms of the fluctuations in the action do 
not disturb the behavior. Since only c and d include the tachyonic field, we estimate the 
orders of their second order terms 5*2 and fourth order terms S4 given respectively as 



^2 = -TpJdP+^x2i(3{c*d~cd*) 

^4 = -Tp f dP+^x\c*d-cd*\^ {A.t 



and show |S'4| <^ 15*21 below. 

TThe appearance of such a shape was already discussed, but was denied in ref.|29| in the case of 
D-strings. 
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For D-p-branes with p > 2, if we include the non-hnear terms in the equation of motion, 
fl2.16|) (with co,fc replaced by its time- dependent typical value) are not the solutions to the 
equation of motion any more. However, since each of /„ and Uk forms a complete set, the 
expansion is still available as a type of decomposition of d.o.f.. So, we substitute ()2.16|) 
for S2 and 5*4 and discuss their amplitudes. Neglecting the contribution of non-tachyonic 
(i.e. non-blowing-up) modes, we have 

^2 = -Tp I dti~2/3)A{tf 

Q f,, ^ fdP-^ ff-^ ff-^., . . 

^4 = -T,J dt^-J (^Co,,,co,.,Co,,3Co,_,,_.,_.3. (4.9) 



since A{t)'^ = J d^^^k/{27ry^^ < |co,feP >• The factor y/^/vr comes from extra normal- 
ization constants of /o's. To carry out the integrations with respect to momenta in ()4.9|) 
is difficult, so we make order estimation of S4, which will later prove to suffice here; first 
we estimate a typical value of A; = Jk^ in condensation of c(x^) using fcav = defined 



at the last part of section 2; it is about 30 ~ 5 percent of yfcmax — monoton- 
ically decreasing during t > lOT^^^ in the case (I) or t > 10T(-^^) in the case (II). Next, 
we choose some appropriate value of fcav and approximate Co,fc. for i = 1,2,3 in ()4.9|) to 
|co,fca^l) defining y4av(t)^ = / (i^^^k/(27r)^^^ < |co^fc^^p >. Then, we can show that the 
behavior of A(t)av is not so different from that of A{t), and we estimate S4 as 

^4 ~ -T,J dt^^^^^{A.At)r. (4.10) 
Thus, the condition for 15*41 <^ 15*21 is 

(27r)P-i/2 



A{t) < 



(4.11) 



We can easily see that for a small angle 6' ^ 2/5 and p > 2, the inequality ()4.11|) is 
sufficiently compatible with ()4.7|1 . The other conditions, especially that of SYM that the 
displacement of branes should be smaller than Ig {A{t) < Is) and D-branes should have 
some low velocity {y <^ 1{= c) at all the value of x), are also shown easily to be satisfied. 
Therefore, the shape of the recombined D-p-brane for p > 2 surely come to have three 
extremes. An example is drown in Fig. 6. This happens essentially due to the localization 
of tachyons around x, since the factor e~^^^^ in ()4.5|1 is directly responsible for the multiple 
extremes. 

For the case of D-strings, the case is a bit more subtle, but we argue that the same 
phenomenon is expected to happen. In this case 5*2 and 5*4 are written as 



^2 = -Toi J rft(-2/?)|go(t)|?ypical 
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Figure 6: The shape of a D-brane with multiple extremes; 2(3 = 0.01 and A{t) = 0.13 



5*4 - -Tdi y"diy ^|co(^)|typical• 
If we substitute Co(t) = A(t) for the condition |S'4| < 15*21, we have 

A{t) < {A7T(3y/\ 
If one consider the equation of motion for co(t), 



2(3co{t) - 2\co{t)\%{t) 



(4.12) 



(4.13) 



(4.14) 



()4.13|1 is replaced by A{t) < {tcPY^^. Anyway, there is, a little though , a range for A{t) to 
be allowed to have three extremes within second order approximation: (7r/9/2)^/^ < A(t) < 
{7ipy/\ This may seem subtle, but if we solve numerically the equation of motion for 
co(t) we can see that |co(t)| really develop to exceed ^critical in ()4.7j) . The other conditions 
of the approximations are again shown to be satisfied. Thus, it is probable that the shape 
of the D-string also comes to have three extremes. 

The physical interpretation of this seemingly queer behavior is as follows; the energy 
released via tachyon condensation pushes the recombined branes away from each other, 
but it is given only to the local part of the branes around x = due to localization of 
tachyon condensation, so, only the part is much accelerated. Though the D-branes have a 
large tension, they also have a large inertia, and when the given energy of the local part is 
large enough, the branes extend, surpassing the tension, to form three extremes. That is, 
localization of tachyon condensation or that of the released energy, and the (large) inertia 
causes the shapes of branes to have three extremes. 

The recombined picture has been obtained, and that the seemingly queer behavior 
can also be explained from physical point of view. Thus, we can say that each diagonal 
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element of Xp+i describes the shape of each D-brane, and that in the gauge where Xp+i 
is diagonahzed (which we denote as "physical gauge"), the correspondence of degrees of 
freedom is direct between open strings and the matrices Xj and A^; for example, (1,1) 
component of corresponds to excitations of open strings whose both ends are on a first 
D-brane, etc.. 

5 Creation of open string pairs 

Let us next discuss creation of open string pairs. We mainly concentrate on case (I) in this 
section: We first discuss the case of D-strings, and then that of D-p-branes with p >2. 

The blow-up of the gauge field appears in the off-diagonal elements of Ap (in this case 
of course p = 1) and a typical value of its electric flux is 



so this excitation should have some physical meaning. What is the physical meaning of 
the blow-up? We argue that it means that open string pairs (each composed of two strings 
with opposite orientations) connecting the two recombined D-branes are created. We will 
explain the reason in the following some paragraphs. 

First, we show that this blow-up correspond to excitations of open strings connecting 
the two branes. Here we would like to confirm some basic facts about SYM as low energy 
effective theory (LEET) of D-branes: As discussed in jnH], each D-brane (Dp^ for 1=1,2) 
has U(l) symmetry (U(l)j-), and open strings starting at Dp^^ and ending at Dpg (which we 
denote as 1— >2 strings) have U(l)-,^ x U(l)2 charges (-1,1), and 2^1 strings have charges 
(1,-1). In the U(2) SYM plus adjoint Higgs theory as LEET of branes, effective fields 
corresponding to the strings also have such charges, and if degrees of freedom (d.o.f.) of 
Dpj correspond to (I,I)-components of the U(2) matrices for 1=1,2 (as in the case of the 
physical gauge), one can represent U(l)^ and U(l)2 transformation as Ui = diag(e*°^, 1) 
and U2 = diag.(l, e'"^) for parameters ai and a2- Suppose we use a unified symbol of the 
U(2) field Am, which is equal to A^ for M = /i (=0, 1, - ■ ■ ,p) and is equal to Xi for M = z 
{=p + 1, ■ ■ ■ , 9). Then, Am transforms under the U(l)-,^ and U(l)2 with constant a/'s as 




This flux has a non- vanishing contribution to the energy density T-Ce\e{x) as 




Am 



( 




) 



( 



ttM 



) 



(5.3) 



■M — 



dM 



CMC 
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so, d.o.f. of Cm correspond to those of 1 — >• 2 strings and d.o.f. of Bm to those of 2 — > 1 
strings. 

In our case it is when Xp+i is diagonahzed (i.e. in the physical gauge) that d.o.f. of 
Dp J correspond to (I,I)-components of the U(2), so we have to see the gauge field in the 
physical gauge. Using Uq in fl4.4|) (and abbreviating Ctypicai to c here), we have the field 
in the gauge A'jyj = UqAmUq'^ + iUodMUQ^ as 



A. 



fix doc*c+c*doc 



c* 




iy2 \c\2 y c 

K = 1 (5.4) 

and the electric flux as 

Fo,M = UoFo,ix,)U,-' = Y ) • ^5.5) 

That is, the flux is invariant under the transformation via Uq] the blow-up arises only in 
its off-diagonal parts and hence, according to the above basic facts, the blow-up of the 
electric flux is due to excitations of 1^2 and 2^1 strings. 

Then, how is the ratio of 1— s>2 and 2—^1 strings? We can see that via charge densities 
Jq^ (for 1=1,2) of current j^^^; the above basic facts means that Jq^^ measures 

the number density of endpoints of open strings incoming to Dp^^ minus that of strings 
outgoing from Dp^ (if it is appropriately normalized), and Jq'^ is equal to —Jq^^ in this 
setting. In the notation of Am in (|5.3|) . Jq^^ is written via Noether's theorem as 

jo\= -j?^) = iiihidoCM - CMdobM) - (&m<9mCo - cm<9a^&o)} 

-2(ao - do)hMCM + (^m - dM){hcM + &mCo). (5.6) 

Substituting the solution (j5.4j) for (j5.6|) . j'q^ is shown to vanish, which means that the 

number density of endpoints of open strings incoming to Dp^ at the point x is equal to 

that of strings outgoing from Dp^ at x. This is also the case with Jq and Dpj. That is, 

open string pairs connecting the two decaying D-p-branes are created, each pair of which 

is composed of 1— i>2 and 2— >-l strings! * 

*We can also look into the configuration of fundamental strings in our solution ()5.4|1 via electric source 
current Jqm{xn) of spacetime NSNS field Bmn by using the results of ref.^Hl, done in ref.pj- Then, 
we can show that Jqm vanishes in this case, which is consistent with j^^^ = Jq^' = and our interpretation 
of it. 
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Furthermore, excitations of the scalar and gauge fields in ()5.4|) also includes other 
components which correspond to such open string pairs. ^ We can extract the components 
by writing all the excitations in terms of U(2) covariant forms in the physical gauge, as 

By T-dualizing with respect to Xp+i direction, DqX'p^^^ and D'^X'^j^^ correspond respec- 
tively to the electric and magnetic fiuxes i^o^p+i and -Fpp+i. Thus, based on the same 
reason as the case of the electric fiux Fq^, we argue that only diagonal elements in ()5.7|) 
correspond to d.o.f. of D-l-branes and the remaining off-diagonal components in ()5.7p 
correspond to d.o.f. of open strings connecting the branes. Consistency check of this 
argument is given in the following way: Non-vanishing part of the total energy density of 
th system is 

^totaKx) = ^Tr{l2x2 + + liD',X'^^,r + (5.8) 

and contribution of the diagonal elements of (j5.7|) to TitotaK^;) is 

Hdiag(a:) = Tdi{2 + 2 ■ ^{d,yf + 2 ■ ^{doyf + 0}. (5.9) 
while that of the off-diagonal elements of (|5.7|) and F^^ is 

Hoff-diag = Tdi{(/3 - 2\dmd\yy' + {Pxf\dod\yy' + \doc\'}. (5.10) 
The first and second term in ()5.9|) are rewritten within second order approximation as 



2^1 + [dxyYi the energy density due to the D-strings' existence (tension), and the third 
term is of course the density of the D-strings' kinetic energy. That is, only 7idiag(a;) just 
gives the energy density that the D-strings themselves have. Thus, it is plausible that the 
remaining energy Tiofr-diag is the energy due to interaction of D-branes, which is consistent 
with our argument that open string pairs connecting the D-branes have the rest of the 
energy. 

We can also estimate roughly the number density of the created open strings; since 
SYM describes the behavior of low energy modes or almost zero modes of open strings, 
the strings connecting the two D-branes should be nearly straight. Since the tension, or 



^1 thank K.Hashimoto for letting me notice these extra contributions. 

"''Further confirmation of our argument that open string pair creation are created between the D-branes 
wiU be given in ref . |54| . 
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the energy density of fundamental strings is constant (=1 in our notation), the number 
density of the open strings (per unit length) n{x, t) is equivalent to the energy density 
divided by the displacement of the two D-p-branes in Xp+i direction. § Thus, it is obtained 
as 



n{x, t) — 



(5.11) 



where we have used dod — \/2pd and Bqc — \/2pc. A distribution of n{x, t) for Qs — IQ ^ 
and A{t) — 0.1 is given in Fig.7, and the situation is roughly depicted in Fig. 8. 




Figure 7: Number density of open strings for A{t) = 0.1 with gg — 10 




Figure 8: A rough picture how open string pairs are created 



§Iii fact, open strings have momenta since the ofF-diagonal elements to momentum density is non- 
vanishing. Here, we ignore the contribution of momenta as an approximation and and estimate the 
number density. 
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In fact, a total number of the fundamental string in Fig. 7 is = / dxn = 50.14, i.e. 
not an even integer. This is due to our partially incomplete analysis, in which we do not 
consider the effect of quantization of charges or fluxes. Actually, n{x, t) is an expectation 
value or a probability of distribution of number density for open string to be created, 
and the real process should happen in a quantum theoretic way; each pair will be created 
rather "suddenly" with the above probability. 

Before moving on to the case of D-p-branes with p > 2, we discuss distribution of the 
released energy. The energy released by the shortening of D-branes is 

-^'released = Tdi j dx 2^1 + {dxy)'^\t=o - Tdi J dx 2^1 + {d^yy\t 

^ Tni J dx{P' - (/3 - 2\d\'f{(3xf/y'}, (5.12) 

and it is used partly to accelerate D-branes and partly to create the open string pairs. 
Then, how much is the ratio of the energy distribution? D-branes' kinetic energy E^i kinetic 
and the energy used to create open string pairs -Eopen are given respectively as 

Edi kinetic = Tdi J dx{doy)'^, -^open = Tdi J dxHos-diag (5.13) 

(Of course it holds -^released = -^Di kinetic + -E'open-) Siuce we kuow all the quantities in the 
formulas of the energies, we can calculate the ratios Edi kinetic /-^'released and -E'open/-£'reieased 
at the time t numerically. The result is 

Edi kinetic/ -^released < 10% 

Eopen/^released > 90% for A{t) < 0.1 

as in Fig. 9. and Fig. 10. 

0.2 
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Figure 9: Edi kinetic /-^released vs. t{xls) 

That is, most part of the released energy is used to create the open string pairs connecting 
the branes] only a little part of it is used to accelerate the D-branes. Thus, from the 
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Figure 10: -Eopcn/-E^released VS. t(xZs) 




viewpoint of D-branes' motion, dissipation of the energy of the branes is very large at 
least at the initial stage of decay. We note that the order of -Eopen in terms of Qs is 
0{{gaY), and asymptotic behavior of -Eopen is -Eopen oc e^v^* or eV^* (its exponent is 
decreasing along with time) for A{t) < 0.1. 

Next we study the case of D-p-branes with p > 2. The main difference is that they 
have extra world- volume coordinates Xp, and hence the condensing tachyon field c (and 
Uq) can have dependence of Xp. Then, in addition to A'^, A'^ and -^p+i in ()5.4|) . Ap 
(/2 = 1, ■ ■ ■ ,p — 1) arises through the transformation via Uq, and both U(2) electric and 
magnetic flux have non-vanishing contribution to the energy. A naive form of Ap is a bit 
complicated, but using the approximation in the footnote under ()4.4|) . we have 



Non-vanishing fluxes and covariant derivatives of X's are 

^*(-''' = ( sic X ) (5-15) 

= ( Jf„ ^^^Y'")- (5-16) 

^ \ (3xdpd/y -dpy J 

in addition to ()5.5|) and ()5.7|) . It should be stressed that the above is U(2) magnetic flux 
but not U(l)'s one on each of the branes. Thus, even if we consider the case of D-3-branes, 
its existence does not mean that D-strings are created. Instead, according to the basic 
facts stated above, these d.o.f. again correspond to the d.o.f. of fundamental strings, and 
jo is again shown to vanish by substituting ()5.4j) and ()5.14j) (with Xp dependence taken 
into account). Therefore, we conclude that in the case of recombination of D-p-branes 
with p > 2, open string pairs connecting the two D-branes are created while they are 
decaying. 
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Let us consider the number density of the string. In the case with p > 2, |9o(ip 
and l^ocp in Tiofr-diag (in ()5.10p ) in the case of D-strings are replaced by \dfxd\ and |9/ic| 
where fi = 0, ■ ■ ■ ,p — 1. But then, if we use << g{k) >> defined in ()3.16|) . it holds 
\dod\^ ^« 2(3 » (/o)V2 and \dfid\^ ~« » (/o)V2, so \df,d\ gives 2(3\d\^ and 
the number density formula ()5.1H) holds again, though in this case it is a density in a 
p-dimensional space. 

The released energy Er{t) is 

Erit) = -Top J d^x [2^1 + {d,yy + {dpyy (5.17) 

Since an extra minus term —{dfiyYlt = — « k"^ » (/o)^MP/l/^ is added to Er(t), the 
density part of the released energy is a bit reduced from that in the formula ()5.12j) . but 
the density of energy used to create open string pairs has the same form as that in (jS.llj) . 
(The density of D-p-branes' kinetic energy is reduced in terms of A{t).) Thus, in the case 
of D-p-branes with p > 2, the ratio of the energy used to create open string pairs to the 
released energy is more than that in the case of D-strings, and we can say again that most 
part of the released energy is used to create the open string pairs. 

6 Creation of closed strings 

Let us assume gg <^ 1. Then, the branes can be regarded as heavy objects, and after 
recombination, they start to move away from each other, though very slowly, but with 
some acceleration in opposite directions. Thus, one can easily expect that the system has 
non-trivial time-evolution behavior of mass quadrupole moments, which causes gravita- 
tional radiation. Below, going out of the framework of SYM, we show that the above 
expectation is true, and evaluate its energy flux, based on the time-evolution of the D- 
branes' shapes. In addition, we discuss creation of closed strings for configurations of 
D-branes beyond the reach of SYM. In this section we mainly discuss the case (I) again, 
but we expect that essential part of the discussion also holds in the case of the case (II). 

Since the coupling of a D-brane to the spacetime metric is GiqM ~ 0{gs) (where 
Gio = 5's/(Mio)^ is the gravitational constant, Mio is a 10-dimensional mass parameter 
and M is the brane's mass), the background spacetime can be regarded as a fiat one 
and the wave we will consider is a deviation around the fiat spacetime. Since we are 
dealing with a motion of D-branes with a relatively slow velocity in order to apply SYM, 
a non-relativistic approximation is applicable. In addition, the motion of the D-branes 
are localized around a; = at least at the initial stage, so the following formulae can 
be applied to this case (see e.g. [H]); the lowest order radiation of gravitational wave is 
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electric quadrupole radiation, and the mass quadrupole moment of a collection of particles 
with masses {m/} in 10 spacetime dimensions is 

lij = ^mA[x[x} - ^S^jix^y] (6.1) 

where 's {i = 1, ■ ■ ■ ,9) represent (space) positions of a particle labeled by I. The angular 
distribution of the energy flux due to the radiation is given by 



'8 

where Ui = (x, — Xi)/\xi — Xi\ where x is the position of the observer and < ■ ■ ■ > denotes 
an average over several wavelength here. We note that in ()6.2|) Xj in lij is to be evaluated 
at the retarded time. To apply the above formulae, we regard the branes as an assembly 
of infinitesimal hyper-planes. Then, lij in the case of branes is translated into 

=Tpj:J ^l + {d^x[nxy^ - l6,,{xin (6.3) 



/=1 



where the index /(= 1, 2) denotes each of the branes. 

Let us first discuss the case of D-strings in the case (I), the simplest case, where D- 
strings moves (decays) toward positive and negative directions of X2- The lij which have 
non-trivial components are 



In = ^Jdx^l + {d,yn8x'~y') 



I22 = dx^l + {d^ynSy' - x') (6.4) 

/33 = /44 = ■ ■ ■ = /99 = / dx^l + {d^yny^ + x^) 

Here we also assume that the branes has no velocity in x direction. If we use 9-dimensional 
polar coordinate system as X2 = cos 9i,xi = sin 9i cos 62, x^ = sin 9i sin 6*2 cos O^, - ■ ■ ,xs = 
sin 9i - ■ ■ sin 67 cos 6^8, xg = sin 6'i ■ ■ ■ sin 9s, we have 



{- sin^ ei{—^y + -(1 - sin^ 01 cos' ^i)(- 



dtd^s Svr '2 ' dt^ ' 2^ dt^ ' 

+ sin^ ^i(cos' 01 cos' 02 - sin' ^2)^^} (6.5) 
dEgj-a,vi _ IOttGio y yd^Ii 2 ,d^l2.2 _ 1 d^Ii d^h^ „n 



where h = / rfxy^l + (9,y) and h = J dx^l + {d^y^xy gs. Since ||^| » \^\ 

(because accelerated most are the part around a; = of the D-branes but dxy{x, t) vanishes 
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at X = 0.), the first terms in ()6.5|) and ()6.6|) are dominant. Using the concrete formula for 
y{x,t) in ()4.5|) . we have the energy flux as 

^ l^l6(20fA{t)' (6.7) 

where A{t) oc eV^* is used. It is of course non-vanishing. Thus, the decaying D-strings 
certainly radiate closed strings as gravitational wave. 

In the cases of D-p-branes with p > 2 in the case (I), only the diagonal elements of 
lii have non-trivial values again due to the integration with respect to x^. To carry out 
the computation of the energy flux of the radiation, it is necessary to regularize the world 
volume of the D-branes, and if we do it, we can easily show that the energy flux is shown 
to be non-vanishing, though the final result depends on the regularization volume. Thus, 
the decaying D-p-branes with p > 2 also radiate gravitational wave. 

In addition, radiation of RR gauge field should also occur after recombination of D- 
branes; each of the D-branes in the initial settings is a BPS brane, so each part of the 
branes should be a source of the flux of RR (p+l)-form gauge field. Thus, RR gauge 
wave should also be radiated if the branes are accelerated and have non- vanishing multi- 
derivatives of their quadrupole moment with respect to time,* which is the case we are 
considering. We expect that the energy flux of the RR gauge wave has about the same 
order as that of the gravitational wave. 

Let us discuss the order of the energy given to the fundamental strings in the case of 
D-strings. In this paper we do not discuss the scale of Mio, and set Miq = Ms{= 1/ls) 
for simplicity. Then, since Gio = and A{tY ~ 0{gs), the flux dEgj.a.^i/dt ~ 0{gl). The 
energy flux of the RR gauge radiation is expected to be the same order as dEgj-avi/dt. 
Taking into account that -Egravi(i^) oc Ait)"^ and the asymptotic behavior of -Eopen, we have 
the ratio of the energy given to the massless closed strings to that of open strings is 

C^^massless closed/dt ^ ^y^^t^^ 9^^^'. (6.8) 
'^-^open/ 

Thus, the energy given to the massless closed string is very small at the initial stage 

compared with that of open strings, but the ratio increases exponentially with a typical 

time-scale T^^) = 1 j ^/l^ 

As for the case (II), the shape of the D-branes in all direction may not be defined due 

to the non-commutativity of the scalar fields, but it is not plausible that cancellation of 

*We note that the lowest order RR gauge radiation is also their electric quadrupole one; the two 
D-branes are regarded to have RR charges of the same signature because if one set = 0, the force due 
to RR gauge field is repulsive. 
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the radiation would happen. So, gravitational radiation and that of RR gauge field are 
strongly expected to occur as creation of closed strings. 

Next, we want to discuss the behavior of the system beyond SYM. If the distance 
between the D-branes becomes larger than Ig, it is difficult to imagine that the created 
open strings extend unlimitedly, because the string has its tension. It also seems unnatural 
if one take into account increase of entropy. Then, what will happen? Since open strings 
are created in pairs it is expected that each pair of open strings connecting the branes 
are cut to pieces to form some closed strings and two open strings, each of the latter of 
which has its both ends on each of the two D-branes, as Fig. 11. That is, such a picture 




Figure 11: 



seems to arise that the decaying (annihilating) D-branes leave many closed strings behind 
(and radiate some of them) while producing many open strings which start from and 
end on each brane. If this is true, the further the two branes go from each other, the 
more possibility there arises to create closed strings. This expectation is consistent with 
tendency of the ratio of dEdosed/dt to dEopen/dt. The dominant mechanism to create 
closed strings in this setting might be the above way. (Please note that if there was a 
single open string connecting the decaying D-branes, the energy of the system was not 
saved by creating a closed string.) 

It is also to be noted that the above picture may be regarded as a generalization of 
Sen's conjecture for DD annihilation [21]; in the case of DD with no topological defects, 
the sum of the tensions of a D-brane and an anti-D-brane are precisely the height of the 
potential of the unstable vacuum and they vanish to be nothing. On the other hand, 
the setting we are considering is a "middle" configuration of DD and parallel two BPS 
branes, the latter of which has definitely positive energy. Thus, our setting also has 
positive energy and hence there should remain something after the decay. According to 
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our speculation, it may be the closed and open strings stated above. 

7 Implications for an inflation model 

In this section we consider implications of the above results of the system for the inflation 
model of the setting of two D-4-branes at one angle in ref.|3]|HI. 

First, we would like to speculate about the efficiency of the reheating of the model. 
According to the analysis in section 5, most part of the released energy via tachyon 
condensation is used to create open string pairs connecting the branes and only a little 
part is used to accelerate the branes, at least at the initial stage of decay. Thus, from the 
viewpoint of motion of the branes, dissipation of the energy is rather large initially, and 
might be large afterward, so the oscillation of the branes around some stable configuration 
might occur only a few times. That is, if this setting is applied to the infiation scenario, 
it is possible that the reheating might be efficient. 

Second, we would like to discuss the mechanism to create "reheated" gauge fields 
and fermions; from the viewpoint of cosmology they should be produced at the end of 
the infiation, but the mechanism how are they produced in this setting of the inflation 
model has not been clarified. We speculate about a scenario of the mechanism based on 
our analysis: Although we are considering a simpler setting model, if a more complete 
model of brane-world is considered, the models should be the one in which at the end of 
infiation, the configuration of branes are put so oriented that they reproduce excitations 
of Standard Model as given in e.g. ref. J3]pni or some others. (A prototype of such 
models are discussed e.g. in ref. p.) In the model, fields hke fermions and gauge fields 
arise as open strings connecting the branes or on each of the branes. Our suggestion is 
that the open string pairs connecting the branes like those in Fig. 8 and the open strings 
on each of the branes like those in Fig. 11 might directly correspond to the "reheated" 
fermions or gauge fields created at the end of the infiation. It would be interesting to find 
a framework describing the process in the region of distance other than SYM or beyond 
SYM, and discuss the above possibilities. 

8 Summary and discussion 

We have investigated what happens in the process after the recombination of two D-p- 
branes at one angle for the two types of concrete initial conditions; case (I):two D-p-branes 
have been parallel until t = but are put intersected at one angle 9 at the instant t = 0, 
and case (II):one of parallel two D-p-branes were rotated by a small angle 9 and are 
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approaching each other very slowly. In summary, the results are the following: 

• Explicit functions of typical amplitudes of condensing tachyonic fields (dependent 
on time and location) have been obtained. 

• Time-evolution of the D-branes' shape after recombination has been obtained, and 
seemingly queer behavior of the D-branes' shape that each brane comes to have 
three extremes has been observed. Physical interpretation of it has also been given. 

• Pair-creation of (many) open strings with opposite orientations connecting the two 
decaying branes has been shown to occur. Expectation value of number density 
of the created open strings has been estimated, and energy distribution into the 
open strings has also been evaluated; it has been shown that most part of the 
energy released via tachyon condensation is used to create the open string pairs, 
and dissipation of the D-branes' energy is large at the initial stage. 

• Creation of closed strings as gravitational radiation from the branes has also been 
shown to occur; its energy flux has been evaluated quantitatively. In addition, a 
strong suggestion has been obtained that closed strings are created out of open 
string pairs after the two branes get further from each other. 

• Speculations on implications of the above results for the inflation model of the D- 
4-branes at one angle have been given: dissipation of the branes' energy (due to 
the creation of the strings) might be large and reheating might be efficient. In 
addition, "reheated" gauge fields and fermion which should appear at the end of the 
infiation, might correspond as the open strings created accompanying the decay of 
the D-branes. 

Among the above, the most important result is the pair-creation of open strings con- 
necting the decaying D-branes, which strongly suggests that closed strings are created 
afterward. These creations in the case of decay of D-p-branes at one angle are in some 
senses like the case of decay of an unstable D-brane (or s-branes|42j) where rolling tachyons 
jlH] are used to discuss production of open string pairs|ll] jlHl and closed strings 

[15] jini ISO] (see also [21] and other references therein) . Compared with the latter 
cases, the analysis used in this paper is only perturbative and can follow the behavior 
of the system for only a rather short time-scale. For further study of recombination of 
D-branes, the approaches analogous to the case of an unstable D-brane may be useful and 
is to be explored. 
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Finally, we discuss the problem of non-commutativity, uncertainty and diagonalization 
procedure (gauge choice) of VEV's of adjoint scalar fields Xp^i and Xg, which appear in 
the solution of case (II) . The commutation relation of the two is 



[X,^^,X,]^z{t)\^''^ Q J. (8.1) 

It is plausible for the relation to give an uncertainty relation 

Axp+iAxg > \z{t)\\dx,t\ oc \zo - vt\A{t)e-'^''\ (8.2) 

That is, this is a solution in which the scale of uncertainty changes depending on location 
and time. 

Let us assume this uncertainty relation. Then, the situation is as follows: If one 
focuses on the positions in Xp+i direction and do not determine those in xg direction, it 
seems possible to discuss the shapes of branes projected on XpXp+i-plane with some high 
accuracy, so, we have argued that the discussion in section 4 holds. As for the "shape" (or 
profile) in xg direction, suppose one diagonalizes Xp^i and describes their positions in the 
direction with the accuracy of the order \d{x = 0)|, i.e. the accuracy one can determine 
whether the branes are combined or not at x — 0. Then, one is able to determine the 
positions in xg direction only with the accuracy of order \z\e~^^^ , which means that we 
cannot tell whether the branes are separated or not at x — 0, and at the point a little 
away from x — 0, one can tell the branes are separated but cannot determine the shape. 
And at the point more away from x — 0, one can tell the shape of the branes projected 
on XpXp+i-plane. 

One may notice that there seems a bit problematic thing: When one takes the gauge 
which diagonalizes Xg instead, no evidence of recombination process can be observed. It 
may be possible to attribute this problem to the gauge choice: On the one hand there is 
a gauge where geometric information of the branes is easy to extract ("physical gauge" 
where Xp+i is diagonalized) . On the other hand, there is a gauge where the behavior of 
fluctuations like tachyon condensation is easy to describe (the gauge Xg is diagonalized). 

In addition, we would like to speculate about the meaning of diagonal elements of 
the scalar field which is not diagonalized; it might have some geometric meaning. In the 
physical gauge, the VEV of Xg is written as 

X, = U,X,U„-^ = ^y=P^ ( ^^-^ ) . (8.3) 

The diagonal elements represent two kink-hke shapes (profiles) along x-direction with 
a width \d\/(3 with both ends approaching ±z{t)/2, which fits our expectation for the 
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branes' shapes projected on XpXg-plane. It may be a merely coincidence, but there might 
be a possibihty that it has some geometric meaning. Let us remember the solution of 
dielectric (spherical) D2/D0 bound states given by Myers in ref.jSSl- concentrate the 
N=2 (two DO) case. The solution is 

^i = A(? J)' ^2 = ^(° o'), Xs = f^(^l \ ^ (8.4) 

where R is the extent of the spherical D-2-brane to which two D-O-branes are bound. If 
one focuses on X3 direction, the two D-O-branes are considered to locate in the north and 
south pole of the sphere in the X3 direction. Then, one cannot determine the positions 
of two DO's in xi and X2 direction, but the zeros in diagonal elements of them may be 
interpreted as center of mass positions or expectation values of their positions. Analogous 
to this case, it might be possible that the diagonal elements of X'g describe center-of mass 
positions of the branes in xg direction though the definite position is uncertain. The above 
point, related with non-commutativity, gauge choice and uncertainty, would be interesting 
and to be explored in the future. 
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Appendix 

A The width of a Gaussian wave function for a mode 
with a time-dependent (negative) frequency-squared 

In this appendix we first review derivation of the propagator (kernel) K of the harmonic 
oscillator with a time-dependent frequency, and then the fact that an initially Gaussian 
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wave function keeps to be Gaussian within WKB approximation, and finally present the 
formula for the time- dependent width of the Gaussian in terms of basis functions. See for 
more detail ref . j37j . 

If one find two independent solutions A(t) and B{t) (which we denote as basis func- 
tions) to the "WKB-approximated equation of motion 

^ = -^{trx{t), (A.i) 

we can immediately write essential part of the propagator [STj in terms of A{t) and B{t) 
up to WKB approximation. The solution representing the classical path of a motion from 
{xo,to) to (a;i,ti) is 

= n A ^ A n [{BiA{t)-AiB{t)}xo + {-BoA{t) + AoB{t)}xi] (A.2) 
= ho{t)xo + hi{t)xi (A.3) 



where we abbreviate A(ti) as Ai, etc. Substituting ()A.3j) for the action S = fl^^ dt m(x^ — 
a;^a;^)/2, and exponentiating the one times i, we have the propagator K{xi,ti; Xo,to) = 
PqlSci "v^^here Sc\ is quadratic in x as 

= IT XiXjaij{ti;tQ) (A.4) 

^ j=0,l 

where 

aijih- to) = ^^[^^ ~ uitfhihj]. (A.5) 

F is a time-dependent normalization but not dependent on Xi, which can be determined 
in several ways. (In this paper we determine the factor F so that the propagator in each 
case K satisfies the Schrodinger equation.) Since the wave function for t = ti is given by 
\l/(xi, ti) = / dxoK{xi,ti; xq, ^0)^(3^0, ^o), the absolute value of an initially Gaussian wave 
function remains to be Gaussian, since the above is only a Gaussian integral. The behavior 
of the system is represented by the time-dependent width of the Gaussian obtained in 
terms of the basis functions as 

Thus, finding the basis functions leads directly to the time-evolution behavior of the 
width and the system up to WKB approximation as far as the the initial wave function 
is Gaussian. We note that normalizations of A{t) and B{t) do not appear in ()A.3jl and 
hence also not in the final result. 
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